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A novel approach for calculating Casimir forces between periodically deformed objects is devel- 
oped. This approach allows, for the first time, a rigorous non-perturbative treatment of the Casimir 
effect for disconnected objects beyond Casimir's original two-plate configuration. The approach 
takes into account the collective nature of fluctuation induced forces, going beyond the commonly 
used pairwise summation of two-body van der Waals forces. As an application of the method, we 
exactly calculate the Casimir force due to scalar field fluctuations between a flat and a rectangular 
corrugated plate. In the latter case, the force is found to be always attractive. 
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Casimir forces between electrically neutral solids are a 
macroscopic consequence of material dependent changes 
in the zero-point vacuum fluctuations of the electromag- 
netic field |l|, [|, [jj. The Casimir effect for two parallel 
metallic plates has been referred to as one of the least 
intuitive consequences of quantum electrodynamics [Q. 
However, this effect can be as well regarded as a macro- 
scopic manifestation of many-body van der Waals forces 
between the particles that form the plates. This point of 
view is supported by Schwinger's approach which consid- 
ers no vacuum fluctuations but only the fields generated 
by the fluctuating dipoles themselves 0. In addition, 
fluctuation induced forces can be observed in a plethora 
of other systems where quantum-, thermal or disorder 
fluctuations are modified by external objects ||. Ex- 
amples of recent interest are as diverse as the thicken- 
ing of Helium films near the superfluid transition || 
forces between vortex matter in anisotropic superconduc- 
tors H , Casimir energy densities in cosmological models 
|| , and interactions between proteins on biological mem- 
branes |J. 

Initiated by the first quantitative verification of the 
electrodynamic Casimir effect by Lamoreaux flfj|| , high 
precision experiments have motivated a resurgence in the 
field of Casimir force measurements in the last few years 
fly]. All experiments confirm the theory for Casimir's 
original flat plate geometry within a few percent of ac- 
curacy. Much less clear is the situation for non-trivial 
geometries. If the separation between the objects is large 
compared to characteristic wavelengths of the material, 
e.g., the plasma wavelength, the force becomes a univer- 
sal function of the geometry with the energy scale set 
by h. However, the entropic origin and collective nature 
of the Casimir force can cause this universal function 
to have a non-trivial and unexpected dependence on the 
shape of the interacting objects. There is little intuition 
even as to whether the interaction is attractive or repul- 
sive as demonstrated strikingly by the positivity of the 
Casimir energy of a conducting spherical shell 0|. In 
particular, the latter result raises the important question 
if repulsive forces can also emerge between disconnected 



objects, being of potential importance for the behavior 
of microelectromechanical systems due to the short-scale 
separations between their mobile parts 13 . Recently, 



non-trivial shape dependencies have been probed experi- 
mentally by tailoring the shape of the plates in Casimir's 
geometry |TJ, [TJ|] . 

Calculations of Casimir forces beyond the simple situa- 
tion of two parallel flat plates are notoriously difficult due 
to the collective nature of the many-body interaction. To 
date, this simple geometry indeed appears to be the only 
case amongst disconnected macroscopic objects for which 
an exact result is available. This is due to the lack of rig- 
orous, non-perturbative methods for calculating the force 
between deformed objects. The simplest and commonly 
used approximation is the pairwise summation (PWS) of 
van der Waals forces S. However, Lifshitz's theory for 
dielectric bodies (with flat surfaces) demonstrates that in 
general the interaction cannot be obtained from such an 
approximation jl6) . In view of the fact that the PWS ap- 
proximation can even give the wrong sign for the Casimir 
interaction, e.g., for a conducting cubic shell fl7fl , it is fair 
to conclude that this approach is to some extent uncon- 
trolled. For corrugated metal plates, the failure of PWS 
for small corrugation lengths has been shown by pertur- 
bation theory with respect to the deformation amplitude 
fll8| . Another perturbative approach is based on a multi- 
ple scattering expansion but has been applied only in the 
limit of large separations between the objects [|9|. How- 
ever, with regard to potential sign changes in the interac- 
tion and experimentally realized separations, it is highly 
desirable to gather information about the perturbatively 
non-accessible regime where the distance between the ob- 
jects is of the same order as the deformation amplitude. 

In this Letter we will introduce an avenue along which 
the shortcomings of the commonly used approximations 
can be bypassed for objects with uniaxial shape modu- 
lations. To illustrate our new approach, we consider the 
force between a rectangular corrugated and a flat plate, 
a geometry which cannot be treated by perturbation the- 
ory due to the edges in the surface profile. Before coming 
to this specific case, let us outline the general approach. 
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In uniaxial geometries the electromagnetic field and thus 
also the Casimir force can be represented as a sum of 
transversal magnetic (TM) and electric (TE) wave con- 
tributions p8| . In order to simplify the analysis, we will 
consider here only TM waves, which are governed by the 
scalar field Euclidean action 



£[*] = i J d 4 A(V$) 2 



(1) 



after a Wick rotation to imaginary time. For TM waves 
we have to impose Dirichlet boundary conditions, $|g = 
on the surfaces S. TE waves can be treated analo- 
gously by imposing Neumann boundary conditions with 
suppressed normal gradient of $. The geometry consists 
of a flat surface in the xy-plane at z = H and a uniaxially 
corrugated surface given by z = h(x) with periodic height 
profile h(x + X) = h(x) of wavelength A. At zero temper- 
ature, the Casimir force per surface area F = —d£/dH 
between the two objects corresponds to the change with 
the objects distance in the ground state energy density 
£ of the field This zero temperature theory is appli- 
cable as long as the de Broglie wavelength of photons, 
At = hc/ksT, is much larger than the mean separation 
between the objects. 

Employing the path integral method to implement the 
boundary conditions at the surfaces 0, the energy 
density is given by £ — —hchiZ/AL with L the Eu- 
clidean system size along imaginary time and 



Z = Z n 



P$ Yl <5[$(r, h(x))]5[$(r, H)]t 



sm/h 



(2) 



with the partition function Z of empty space with no 
objects, and r = (ict,x,y) = (X ,Xi 7 X 2 )- The delta 
functions can be rewritten as a functional integral over 
an auxiliary field on each surface so that after integrating 
out the new action is quadratic with matrix kernel 



M(r, r') 



G[r - r', h{x)-h{x')) G[r - r', H-h(x')] 
G[r-r',H-h(x)] G[r-r',0] 



(3) 



and Green's function G(r, z) = (r 2 + z 2 )- 1 /4vr 2 §|. The 
Casimir force per surface area A between the objects is 
then given by 

tic nc 
F = -^d H lndet M = -— Tr (M^dnM) . (4) 
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This expression is suited for a non-perturbative approach 
since the force must be finite, in contrast to the energy 
density £ which contains without an explicit cutoff also 
infinite but H independent contributions from the indi- 
vidual surfaces. 

In the more conventional approaches, the zero-point 
energy is calculated in terms of the allowed frequencies 
u> n of the space between the objects, using £ — \ J2 n ^ UjJn 



and by applying subtraction schemes in order to obtain 
a regularized energy which in turn yields the force. An 
advantage of our approach is that the frequencies need 
not to be calculated explicitly but the force is directly ob- 
tained in terms of the free space Green's function without 
any necessity for regularization schemes. Of course, the 
functional inverse of the matrix kernel M(r, r') in Eq. 
(Q) is in general difficult to evaluate. In the following, 
we will demonstrate how this problem can be tackled for 
uniaxial geometries. Due to the symmetry of the geome- 
try, it is useful to consider the matrix kernel M(p,q) in 
momentum space. The periodicity of the surface profile 
along the X\ direction and uniformity in the XqX 2 plane 
impose the general form 

oo 

M(p,q) = N m (p ± , Pl )(2n) 3 S(p ± + q ± ) 

TTl — — OO 

x 5{pi +qi + 2-n-m/X) (5) 

on the kernel. This Fourier decomposition defines the 
matrices N m (p±,pi) with pj_ = {p Ql p 2 ) and p± = |p_ij. 
Thus the matrix M(p, q) has non- vanishing entries only 
along the diagonal and its periodically shifted analogues. 
An important observation is that matrices of this struc- 
ture can be transformed to block-diagonal form by ap- 
plying row and column permutations. Physically, the ex- 
istence of such a transformation is due to the non-mixing 
property of the kernel M(p, q) for Fourier modes whose 
moments differ by non- integer multiples of 2ir/\. For the 
time being, let us assume that the objects have the exten- 
sion W along the X\ direction, leading to a discrete set 
of momenta p\. Then the transformed matrix Af(p, q) 
consists of block matrices Mj along the diagonal with 
j = 1, . . . , N = W/\ — 1. Since this matrix form can be 
always realized by an even number of permutations, we 
get 



detM 



II 



det Mi 



(6) 



The block matrices Mj can be expressed in terms of the 
matrices N m of the decomposition in Eq. (||). For dis- 
crete momenta in x direction, the matrices Mj can be 
parametrized by integer indices fc, I = — oo, . . . , oo, lead- 
ing to the result 



M. 



,«(px,q±) - (27r) 2 S(p ± + q ± )B kl ( P± , 2nj/W) (7) 



itli the matrix Bki given by 



B kl { P± , Pl ) = N k _i{p ± , Pl + 2ttZ/A). 



(8) 



Moreover, from Eq. (0) we obtain easily 8h (In det M) = 
^2j = i Tt{Mj diiMj). Thus the separation dependent 
part of the ground state energy £ is given by the sum of 
the individual energies of the decoupled subsystems which 
are described by the matrices Mj. Next, it is useful to 
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define the function 

g(px,Pi) = tr (B _1 (p± > pi)aHS(p±,pi)) , (9) 

where the lower case symbol tr denotes the partial trace 
with respect to the discrete indices k, I, cf. Eq. at 
fixed p±_ , p± . It follows from the structure of the matrix 
Bki that g(p±,p\ + 2irn/\) = g(p±,Pi). In addition, 
choosing the height profile to be symmetric, h(—x) — 
h(x), imposes the condition N m (p±, —p±) — N_ m (p±,pi) 
on the matrices N m which in turn leads after row and 
column permutations of the matrix Bki to g(p±, — Pi) = 
g(px,Pi)' Taking at this point the thermodynamic limit 
W, N — > oo, transforms the sum over the energies of the 
N subsystems into an integral. Thus we obtain, after 
taking the total trace over all degrees of freedom, the 
final result for the Casimir force density, 



F = - 



he 

4^ 



dp±pj 



r/A 



dpig(p±,Pi), (10) 



where we have used the symmetries of g(p±,pi). For 
a given Fourier decomposition of the kernel M into the 
matrices N m , this formula together with Eqs. (@), (g) 
yields the exact result for the force between the objects. 
This result has the advantage of being particularly suited 
for a numerical analysis as will be shown below. 

In the following, we apply the general formula in Eq. 
( |l0| ) to obtain a non-perturbative result for the force be- 
tween a flat and a rectangular corrugated surface. We 
choose the height profile on a basic interval as h(x) = a 
for \x\ < A/4 and h(x) = -a for A/4 < \x\ < A/2 so 
that H is the mean distance between the surfaces, see 
Fig. |l|. This profile allows for an analytical decomposi- 



H 



2a 



FIG. 1: Cross section along the zz-plane of the two-plate 
configuration with translational invariance in y direction. 

tion of the kernel into matrices N m (p±,pi); details of the 
calculation will be published elsewhere. The result is 



No = 2p 



$ o(p) ^-cosh(ap) 



2p 



cosh(ap) 




_i_ 

2p 



2p 



■ sinh(ap) 



V 2p m smh(ap m ) 



for odd m ^ 0, and 



N m = 



*m(p) 




(11) 



for even m ^ with p m = \jp\ + (pi + 27tto/A) 2 , c m 
2(_l)(IH-i)/a/ 7r | m | an d 

_ (-1)W2 ~ f\p ± , Pl + (2n-l)2n/X] 
*rn[P)- ^ 2 2^ (2n-l)(m-2n+l) ' 



n— — oo 



where f(p±,Pi) — (exp(— 2ap) — I) /p. In order to calcu- 
late the function g(p±,pi), we truncate the matrix Bki 
symmetrically around (fc, I) — (0, 0) at order M so that k, 
1 = -(M-l)/2,... ,(M-l)/2. This defines via Eq. (g) 
a series of approximations gAi(p±,Pi) which converges to 
g(p_L,Pl) for M -> oo. 

Before using the general form of the decomposition, let 
us consider two instructive limiting cases. In the limit 
A/a — > 0, one would expect that <7m(p_l,Pi) converges 
very rapidly since the contributions from matrices N m 
decrease with m, cf. Eq. (||). To test this guess, we make 
use of the fact that the matrices N m (px, Pi +2%l/X) sim- 
plify considerably for A/a — > so that the series gu can 
be explicitly obtained order by order from the truncated 
matrix Bki- Using Eq. (0), we get 



5m(p_l,Pi) 



2p(l+e- 



!) 



1 + e -2ap_ 2e 2(ff- 

2p 



for M = 1 
for M > 3 



(12) 



for positive a. This result indeed shows a rather rapid 
convergence; the series gu is even invariant for M > 3 
with increasing dimension of the matrix Bki ■ Using Eq. 
(|l0|), we get the force density 



he 



480 (H - a) 4 



(13) 



Physically, this result appears to be quite natural since 
the most contributing modes of the fluctuating field can- 
not probe the valleys of the corrugated surface in the 
limit A/a — > 0. Thus, the surfaces experience a force 
which is given by the force between two flat surfaces at 
reduced distance H — a, corresponding to Eq. (13) [5j| . 
In the following we call this limit the reduced distance 
(RD) regime. For small a/H, the correction ~ hc\a\/H b 
to the flat surface result is non- analytic in a, and thus 
cannot be obtained in perturbation theory 

Next, we consider the limit X/H — + oo. In the limit of 
small surface curvature, the pairwise summation (PWS) 
of renormalized van der Waals forces is justified p3[ . For 
a general surface profile, such a summation corresponds 
to a spatial average over x of the Casimir energy for two 
flat surfaces, taken at the local distance H—h(x) fisj . For 
the geometry of Fig. [l], this procedure obviously yields 



Foe — 



7r 2 he 



1 



1 



480 2 \{H + a) 4 (H - a) 4 



(14) 



since 50% of the local surface's distances are H+a, H — a. 
each. 
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FIG. 2: Relative increase of the Casimir force density F 
for the geometry of Fig. Pll compared to the result Ffl at = 
— 7r 2 /480 x hc/H 4 for two flat surfaces. 
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FIG. 3: Crossover between the PWS (F -> and the RD 
(F — > .Fo) regimes for fixed H/a. Inset: Scaling behavior of 
F for the approach to Fo for A — > and to Foo for A — ► oo. 
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In order to see how the previously discussed limits fit 
into the complete theory as described by Eq. (|lO|), one 
has to resort to a numerical analysis. The recipe for this 
analysis is straightforward: At fixed order M, the trun- 
cated matrix Bm is calculated from Eqs. (|ll|), then its 
inverse is constructed to obtain the M th order approx- 
imation qm to the function g. Finally, numerical inte- 
gration in Eq. (|lo| ) yields a corresponding series of force 
densities which can be extrapolated to the final result for 
Al — > oo. The results are summarized in Figures |^ and 
^. For fixed A/a, the relative increase S — F/Ffl at — 1 
of the force compared to the force between two flat sur- 
faces is shown in Fig. |[ At small A/a the exact result 
Fo of Eq. Jl3| ) is recovered. For larger A/a there is a 
crossover between two scaling regimes: For H ^> A the 
RD scaling 5 ~ H^ 1 remains valid with an amplitude 
which decreases with increasing A/a. For H <C A the 
PWS regime is entered with 5 ~ H~ 2 . As a consequence, 
the force is always attractive with F$ as upper and F^ 
as lower bound. Fig. || shows the crossover between the 
PWS and the RD regimes for fixed separations H/a. In- 
dependent of H/a, the crossover appears at A/a m 10. 
The two limits are approached as (Fo — F)/Ffl at ~ A/a 
for A — > and (F — F 00 )/Ffl at ~ a/A for A — > oo, cf. the 
inset of Fig. |[ 

The reported results may stimulate new developments 
in a large class of systems with fluctuation induced in- 
teractions. Most intriguing might be the possibility of 
a repulsive force for the full electromagnetic field due to 
anomalous contributions from the TE modes fl8|| . Fur- 
ther extensions may include many-body interactions for 
more than two objects and the dynamic Casimir effect. 
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